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Abstract 

We present a detailed thermal and structural analysis of a 2m class solar telescope mirror which is subjected to a 
varying heat load at an observatory site. A 3-dimensional heat transfer model of the mirror takes into account the 
heating caused by a smooth and gradual increase of the solar flux during the day-time observations and cooling 
resulting from the exponentially decaying ambient temperature at night. The thermal and structural response of two 
competing materials for optical telescopes, namely Silicon Carbide -best known for excellent heat conductivity and 
Zerodur -preferred for its extremely low coefficient of thermal expansion, is investigated in detail. The insight gained 
from these simulations will provide a valuable input for devising an efficient and stable thermal control system for the 
primary mirror. 
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1. Introduction 



At the beginning of 17th century, Galileo Galilei turned his small telescope towards the sun and noticed the dark 
regions on sun's surface now called sunspots. The larger telescopes built subsequently revealed that sunspots have dif- 
ferent structures like dark region called umbra and less dark with filamentary structures called penumbra. As the tech- 
nology improved in 20th century, bigger telescopes were built to observe the characteristic f eatures which include th e 
light bridge across sunspots ( VazquezUl973), moving magnetic feat ures around sunspots (IHarvev & Harvevi 1 19731) . 
supergranulation network and its elements (ISimon & Leighton[|l964l) etc. Most of these observations have been done 
using about 50 cm class telescopes or smaller. From the theoretical calculations, the size of the elementary building 
blocks (size of the flux tube) of magnetic structures, for example, is estimated to be about 70 km. A large aperture 
telescope with diameter of 2m or more is therefore, necessary to resolve important and intriguing features of the solar 
atmosphere at smaller scales. Additionally, the photon starved polarization measurements, crucial to estimate the 
magnetic fields distribution and strength in the active regions on the sun, are hard to come by from small telescopes. 

The next generation, large aperture (2m and above) solar telescopes equipped with adaptive optics and supported 
by advanced backend instruments, will be able to probe the solar atmosphere at unprecedented details and thus help 
solving many outstanding problems in solar physics. Some of the notable solar facilities planned for the future are : 
a 4m advanced technology solar telescope (ATST) to be build by National Solar Observatory (Ke il et all 12009), 
Natio nal Large Solar Telescope (NLST) proposed by Indian Institute of Astrophysics, INDIA dSingh , 2008; Hasanl 
20 lob, and a soon to be com missioned 1.5 m Gregor telescope by a consortia led Kiepenheuer-Institute, Germany 
(IPeter & von der Luehel |2006 ) . 

The design and construction of a solar telescope is considerably different from a typical night-time telescope. 
Among others, the most formidable task is to build an efficient temperature control system capable of handling the 
large amount of heat generated by the intense solar radiation. The immense difficulty of the heat management problem 
has been a major bottleneck in constructing large aperture primary mirror for solar studies. Evidently, by the time 
largest solar telescope (4m ATST) starts its operation, the night-time astronomy would have leap-frogged to next 
generation of 20-40m class telescopes. 
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Besides the direct heating caused by light absorption, the diurnal temperature variations of ambient air are also 
crucial in driving the thermal response of the telescope mirror. The rise or fall in mirror temperature with respect to 
ambient, adversely affect the imaging performance of the telescope in two major ways. First, when the mirror is at 
different temperature than the surrounding, a thin layer with a temperature gradient is formed closed to the mirror 
surface and air. The temperature fluct uations in the surface layer lead to the variations in the refractive index of th e 
air that produce wavefront aberrations ( Lownel 1979 Salmon, 1992 : Barr & Livingston , 1993 ; Wood & Ryan , 1995 ). 
The final image is blurred as the telescope beam passes through this layer twice. The second detrimental effects of 
heating or cooling is the structural deformations of the mirror arising from thermally induced st r ess and temperatur e 
inhomogeneities within the substrate volume ( Pearson & S tepp. 11987k ISpitzer & Bolevl 1 19671 iBurge et all [l999). 
Any deviation from the nominal temperature distorts the mirror surface from its desired shape. Thermally induced 
surface imperfections therefore, deteriorate the final image quality of the telescope. 

To minimize the effect of mirror seeing, the temperature difference between the ambient air and the mirror surface 



should be maintained within a very narrow (< ±2°C) range (lEmde et all 2004; Dalrvmple et al„ 2004). At the same 



time, the temperature difference across the reflecting face of the mirror should not exceed +0.5°C. For night-time 
astronomy, the t emperature of the mirr or and telescope dome can be effectively controlled by air conditioning system 



during the day (Mi vashita et all 12003). Blowing cool air over the mirror surface is helpful in mitigating the mirror 
seeing effects. Forced air flow breaks the thermal plumes closed to the mirror surface and thus homoge nize the 
tempe rature gradients which otherwise drive the refractive index fluctuations leading to image degradation (ILowne , 



19791) . The induced air flow also improves the convective heat transfer rate and thereby reducing the temperature 



difference between the mirror material and the surrounding air. For effective cooling, the air temperature and speed 
have to be regulated according to the temperature profile of the mirror surface. However, the accurate measurements 
of the temperature of the mirror require an array of temperature sensors placed close to the reflecting face of the mirror 
which is not just difficult but impractical during the scientific observations. In such a case, numerical simulations are 
the only viable solution to predict the temperature profile of the mirror surface. 

In another approach, a 'cool reservoir' is created by lowering the the mirror temperature well below the temper- 
ature expected for the night observations. A near thermal equilibrium with ambient is accomplis hed in a resistive 
heating by passing electrical c urrent through reflective coating on the front surface of the mirror ( Bohannan et al.l 
2000; Greenhalgh et all 1 1 9941) . This approach does not quite suit the solar observations as the primary mirror is di- 



rectly heated by sun's radiation. The effect of mirror seeing is even more severe if the same telescope is later switched 
to night-time observations. To avoid large heat loads and thermal gradients, ideally a substrate material should have 
low heat capacity and high thermal conductivity. 

Among traditional materials used in primary mirror substrates, ultra low expansion ( ULE) glass ceramics s uch as 
zerodur and fused silica have gained wide acceptance among astronomical community (IDohring et all 120071). The 
semiconducting SiCis another important ma t erial that is beginning to find a niche in astronomy (IBender & Ewing , 
19931: iRiti et aUll99dlPetrovskv et al.Ul994tlTsuno eTafl 120051: Matson et all 12008). In addition to its low thermal 



expansion coefficient, high hardness and rigidity, SiC has exceptionally high thermal conductivity which is about 
two order of magnitude higher than the other materials. For large size telescopes, the lightweighted mirror geometry 
is preferred to reduce the overall weight and thermal inertia of the system. The reduced mass facilitates efficient 
and faster cooling necessary for the mirror to reach the thermal equilibrium with surroundings quickly. However, 
lightweighting is a complex, time consuming and often an expensive process. With a good thermal control system, a 
classical primary mirror of diameter not exceeding 2m could possibly be used without lightweighting. 

For a given shape and geometry of the mirror, the analytical solutions become intractable. Therefore, in tradi- 
tional engineering practice, the thermal respon se of such a mirror blank is studied by finite element methods (FEM) 



(ICho & Hansen . 1995; Dalrvmple et all 120041) . The heat transfer problem described by a diffusion-type partial dif- 



ferential equation, is numerically solved to obtain a steady-state solutions for a constant temperature typically found 
at telescope location. In addition, the mirror heating is modeled assuming fixed heat flux and fixed ambient air tem- 
perature which serves as the boundary conditions. The time scale for a mirror to reach a thermal equilibrium is of 



the or der of few days while the ambient air temperature changes significantly in course of the day (IBohannan et al 



2000). This means, the usual mirror substrate with large heat capacity can never reach the thermal equilibrium with 



ambient unless some external cooling is used. The static analysis is not enough to determine the level of accuracy and 
other parameters of interest (e.g. level of wind flushing or ventilation points etc) necessary for designing an efficient 
thermal control system. A complete time dependent solution of 3D heat transfer problem under the varying solar flux 
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and ambient temperature is therefore necessary to evaluate the mirror performance, a problem which we address in 
the present paper. 

In this paper, we have solved the 3D heat transfer and structural model that takes into account the ambient thermal 
conditions existing at the telescope site. The location dependent solar flux and ambient heating model is incorporated 
into commercially available finite element analysis software to investigate the time dependent thermal and structural 
response of 2m class primary mirror. These numerical simulations are carried out at four different ambient temperature 
range, representative of mild to extreme thermal variations that may occur at different observatory locations around 
the world. 



1.1. Formulation of heat transfer problem 

The general heat flow problem involves all three modes of heat transfer, namely conduction, convection, and 
radiation. The reflecting face coated with a thin metallic layer absorbs about 10-15% of the total solar flux incident 
on the mirror surface. The heat generated as a result of absorption of optical energy is diffused to other parts of the 
mirror via conduction and partly lost to the surroundings via convective and radiative processes. The temperature of 
the mirror rises when the heat produced by light absorption exceeds the overall he at loss to surroundings. Th e general 
heat transfer problem can be cast into the partial differential equation of the form (llncropera & DeW itt. 2001 ): 



P C p ^ - V • (*• VT) +p Cp uV • T = Q 



(1) 



where p, T, k and Cp are density, temperature, thermal conductivity and heat capacity of the material. The terms 
u and Q represent the velocity field and heat source, respectively. Most of the material properties are temperature 
dependent. The boundary condition for the heat flux to be maintained at the surface can be written as: 



n • (*Vr) = q + h(T amh -T) + ecr(T 4 



T 4 ) 



(2) 



where qo is the inward heat flux at the mirror surface, r arm -, is the ambient air temperature, T s is the effective 
radiation temperature of the surroundings, e is surface emissivity, <x is the Stefan-Boltzmann constant and h is the heat 
transfer coefficient indicating the rate at which heat is exchanged between the mirror and the surroundings. The 2nd 
and the 3rd terms in Eq.(2) are statements of Newton's law of cooling and Stefan-Boltzmann's law, accounting for the 
heat loss by free convection and radiation, respectively. 

1.2. Modeling the diurnal temperature cycle 

The diurnal and annual temperature variations is one of the important criteria for selecting a suitable observatory 
site. It also plays a vital role in determining the degrading effects of atmospheric turbulence, mirror seeing and 
structural stability of the telescope. The earth's surface is mainly heated by the irradiated solar energy during the day . 



The sun's light also contains the spectral signatures of both solar and terrestrial atmospheres (Walla ce et al 1994 1. 



The amount and duration of the solar radiation 7(Wm 2 ) reaching the earth's surface is given bv dlaballll983ll 

I = Iq cos(Z) 

where Iq is the flux amplitude when sun is at zenith, and Z is the solar zenith angle which is given by 

Z = arccos(sin 6 sin <p + cos 6 cos cp cos H) 



(3) 



(4) 



where <p is latitude of the place, 6 is declination angle of the sun, H - (tt/12) t is the solar hour angle with respect to 
noon and t is the local time. 

The exact nature of daily temperature variations at a given location, however, depends on the local weather condi- 
tions (pressure, wind speed, h umidity etc), surface topography, soil type, vegetation, and the presence of water body 



etc. |Gottsche & Olesen (2001) proposed a simple physics-based model to describe the thermal heating of the Earth's 
surface in cloud free conditions. According to this model, the day-time rise in temperature T\{t) due to the sun and 
the night-time cooling Ti{t) can be described by 



T\(i) = sTq + s T a cos 
T 2 (t) = s(T + ST) + s 



-(t - r m ) 



t < t. 



-(*-T m ) 



-6T)e 



for 
for t>t 



(5) 
(6) 



For heating part, the choice of harmonic term in Eq.(5) is based on the solution of equation of thermal diffusion 
(ICarslaw & Jaegenll959l') . while the exponential term in Eq.(6) accounts for the temperature fall in accordance to the 
Newton's law of cooling. The constant t„, is included to account for the thermal time lag between the peak solar flux 
and the peak ambient temperature during the day. The meaning and values of various parameters in Eq.(5) and Eq.(6) 
are listed in Table Q] These values are arbitrarily chosen to simulate a typical temperature profile for - 10°C ambient 



Table 1 : Model parameters and the values for Eqs[6] 


Parameter 


Meaning 


Values 




temperature amplitude 


6.8 °C 


T 


residual temperature 


3.4 °C 


ST 


T - T(t -> co) 


-3.5 °C 


to 


half period oscillation of cosine term 


12 hrs 


K 


attenuation constant 


3.5 hrs 


T m 


time when temperature maxima is reached 


14:00 hrs 


u 


time when temperature attenuation begins 


18:00 hrs 


s 


scaling factor 


1.0 



temperature variations. However, a location specific parameters can be obtained by model fitting the temperature data 
from the observatory site. The Eq.(3)-Eq.(6) prescribe the necessary boundary conditions and initial conditions in our 
heat transfer model. The FEM model, thus, has all the necessary input for a realistic thermal analysis of the mirror. 



1.3. The FEM simulations 

For these studies we have chosen a Gregorian type primary mirror design which is most suitable for the solar 
observations. The Gre gorian configuration allows the heat stop to be placed at the prime focus which is located before 
the secondary mirror ( von der Liihel 2009 1. This position is convenient to effectively block the unwanted radiation 
and prevent the subsequent heating of the secondary mirror and other optical components in the beam path. We chose 
a reflecting primary mirror with parabolic shape defined by z = r 2 /4/. The mirror diameter (2r) is 2m, focal length 
/ is 4m and radius of the central hole is 0.21m. The edge-on geometry of the soli d mirror is shown in Fig. [1] This 
configuration is similar to the one under consideration for a 2m class NLST mirror dHasanl l2010HSingh[ bOOsl The 
mirror is primarily designed for high resolution solar observations, but may partly be adopted for some limited but 
valuable scientific observations at night. 



0.25 m 




0.5m 



Figure 1: Mirror geometry for heat transfer model. Three reference points are shown at radial distance r = 0.5 m and different axial positions z 
along the mirror thickness. 

We used COMSOL multiphysic^] as our FEM tool to study the thermal and structural characteristics of the solid 
mirrors. A free mesh containing over 34000 tetrahedral elements was created on the 3D geometry of the mirror. The 
reflecting surface of the mirror was modeled with a 100 nm thick aluminium coating layer. The thin aluminium layer 
requires a much finer mesh which usually take longer computation time. Therefore, we have used the highly conduc- 
tive layer approximation which does not require the resolved mesh for thin aluminium layer. This approximation is 



1 COMSOL Multiphysics Pvt. Ltd., http : //www . comsol . com/ 
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valid since the thermal conductivity of aluminum is higher than the substrate material (see Table [2] for comparison) 
and temperature across the aluminium layer thickness is nearly constant. In other words, any temperature gradient in 
the / = 100 nm thick aluminum layer gets balanced over a minuscule time lag t ~ l 2 /(k/pCp). 

As a generic case, the solid telescope mirror in the FEM model is exposed to solar radiation continuously between 
the sunrise (6:00 hrs) and sunset (18:00 hrs). Light rays are assumed to fall normally on the reflecting surface of 
the mirror. The peak solar flux reaching the earth surface at noon is assumed to be 1000 WrrT 2 . The thermal time 
lag between the peak solar flux and ambient air temperature is taken to be about 2 hours. Which means T am <^ 
continues to rise well past local solar noon, until 14:00 hrs. Eq. (3) is the source of inward heat flux generated at 
the front surface of the mirror during day-time heating t < t s . We have considered free convective heat loss from 
the mirror boundaries (front, back and sides) having heat transfer coeffici ent h = 5Wm~ 2 °C~'. The re l evant thermal 
and mechanical p roperties of the mirror materials are listed in Table l2l dLienhard & Lienhard , 2008 ; Pierre [2002 



Ahmad. A.L 119971) . The FEM model also includes their temperature dependence. 



Table 2: Physical properties of materials at the room temperature. 



Property 




Value 




Units 




SiC 


Zerodur 


Al 




Thermal conductivity (k) 


118 


1.46 


237 


wr'm- 1 


Specific heat capacity (Cp) 


672 


821 


910 


Jkg^KT 1 


CTE (a) 


3.7 


-0.08 


25 


10~ 6 K l 


Material density (p) 


3200 


2530 


2700 


kgirT 3 


Young's modulus (E) 


414 


90.3 


70 


10 9 Pa 


Poisson's ratio (v) 


0.14 


0.24 


0.33 




Emissivity (e) 


0.9 


0.9 


0.05 





The finite element structural model was used to predict the temperature induced mechanical deformations of 
the mirror blank. The mirror substrate will expand or contract during the heating and cooling cycles. Thermo-elastic 
response of the mirror can be examined at any desired instant by mapping the corresponding temperature patterns from 
the thermal model to structural model. The structural FEM model solves for material deformation for 3D geometry 
of the mirror. A fixed boundary condition was applied to the bottom face of the mirror to avoid the rigid body motion. 
The surface nodal displacement was computed for the remaining boundaries that were left free of external constrains. 
No smoothing, transformation or interpolation of the thermal data was necessary as both thermal and structural models 
shared a common mesh on COMSOL software. Simulations were also performed at other ambient temperature range, 
0-5°C, 0-20°C and -10-0°C which correspond to the scaling parameter s = 0.5, s - 2 and s = — 1, respectively, in 
Eq.(5) and Eq.(6). 



2. Results and discussions 

The 3D transient heat transfer model for Zerodur and SiC mirror was solved in two separate stages. The first stage 
run computes the solution for the day-time (6 : 00 < t < 18 : 00) heating and second stage run gives the solution 
for the night-time (18 : 00 < t < 6 : 00) mirror cooling. For each incremental time step (Ar=100s), the ambient 
temperature input was computed from Eq.(5) and Eq.(6). Shown in Fig. 2(a) is the temperature evolution of three 
reference points located on (I) the front surface of the mirror (r = 0.5m, z = 0.203m), (II) the mid-plane of the mirror 
(r = 0.5m, z = 0.101m) and (III) the bottom of the mirror substrate (r = 0.5m, z = 0.0m). For comparison, the model 
input for the solar flux and 0- 10°C ambient temperature curves are also drawn in Fig. 2(a). The grey dotted vertical 
line demarcates the heating (t<t s ) and cooling (t>t s ) regions. 

2.1. Code initialization 

Initially at t = 0, the mirror and the ambient temperature was assumed to be 0°C. Variations in ambient temperature 
and Zerodur mirror for the first 24 hrs simulation run are shown in Fig. 2(a). As the day progresses, the different parts 
of the mirror substrate exhibit dissimilar thermal lags in response to varying solar flux and ambient temperature. Also, 
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for a significant part of the day, no part of the mirror exceeds the temperature of the ambient air. The hot ambient 
also contribute to the mirror heating besides the apparent heating caused by the light induced absorption at the front 
surface. This is evident from the reference point at the mirror bottom (z = 0.0m) which remains at higher temperature 
w.r.t. the midplane point (z = 0.101m) until late afternoon hours when the temperature of the mirror interior is raised 
by the heat conduction from the front surface. After the sunset, the ambient air cools at a relatively faster rate. 




i 18 21 24 3 20 40 60 80 100 120 140 160 180 

Time (hrs) Time (hrs) 



Figure 2: Variations in ambient temperature and Zerodur mirror with time, (a) Simulation run for the first day. The day-time Solar flux curve 
is included to show the thermal lag between mirror temperature and ambient temperature. The scale for the solar flux is shown to the right, (b) 
Simulation ran for the seven days observing cycle. 



The slow fall in mirror temperature due to high thermal inertia (the mass M and total heat capacity MC P of the 
Zerodur and SiC mirror being ~ 1672 kg and ~ 21 15 kg, and 1.37x 10 6 JK 1 and 1.42x 10 6 JK 1 respectively), cannot 
keep up with the rapid fall in ambient air temperature. The mirror and ambient temperatures differs by about 2-4°C 
towards the end of first day simulation run. This has direct implications in the subsequent simulation for the next day 
observation cycle. The initial temperature conditions for the mirror in the next-day run would now be considerably 
different (2-4°C higher than surroundings) than what they were for the start of the first day of observation. Fig. 2(b) 
shows the simulation run for seven days cycle. For clarity, the temperature curves for the front surface reference 
point (z = 0.203m) and the ambient r am ^, alone are shown. Since the mirror does not cool to the surrounding 
temperature during night, there is an overall temperature increase of about 1°C above the ambient. Except for the 
first-day simulation run, the thermal response of the mirror is identical for all other days. Mirror however, never 
attains thermal equilibrium with the ambient. From computing view point, the first-day simulation is more like a code 
initialization phase meant to produce a definitive starting temperature within the mirror for subsequent simulations. 
The code initialization is implicit in all subsequent results presented for Zerodur and SiC in remaining part of the 
paper. 

2.2. Nonuniform thermal response 

Temperature of the mirror rises during the day-time observations when the heat produced at the front surface is 
conducted to the mirror interiors. The mirror remains significantly warmer during the night as the heat is slowly 
diffused now from mirror interiors to boundaries. On an average, the heating rate of the mirror is about 200W during 
the day-time. Considering the 10% absorption, the total heat deposited in the mirror during every 12 hours run is 
~ 8.6 x 10 6 J. As evident from the sliced plane visualization of the mirror volume shown in Fig. 3, the heat distribution 
within the Zerodur mirror is highly nonuniform. The low thermal conductivities of Zerodur is primarily responsible 
for creating axial and radial temperature nonuniformities inside the mirror volume. The shifting temperature gradients 
can be inferred from the fact that a predominant heat concentration is observed in the upper half of the mirror during 
the day as in Fig. 3(a), while the overnight cooling effectively lowers the temperature leaving the core a little hotter 
than the rest of the mirror. But from the seeing perspective, it is the temperature difference (hence forth denoted by 
AT) between the mirror and the ambient air which is of a greater concern. The computed variations in AT for three 
reference points in Zerodur and SiC mirrors are shown in Fig. 4. For Zerodur, the temperature difference between 
the reflecting surface and ambient could change from 1-5°C within 24 hours but it is always above the ambient. 
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Min: 1.44 °C 



Figure 3: Nonuniform temperature distribution inside the Zerodur mirror (a) due to solar heating during the day and (b) convective and radiative 
cooling during the night. 

Therefore, to minimize the seeing effects, the front surface of the mirror has to be cooled continuously during the 
telescope operation. The middle and bottom section of the mirror can have temperature less than ambient during the 
day. The thermal lag between different parts of Zerodur material is because of its low thermal conductivity. 




Time (hrs) Time (hrs) 

Figure 4: Differential temperature variations evaluated for a reference point on top, middle and bottom of (a) Zerodur (b) SiC mirror. The ambient 
temperature range is - 10°C. 

In comparison to Zerodur, SiC has high thermal conductivity. Any heat injected into the mirror is quickly diffused 
into the entire substrate. As a result the entire SiC mirror block attains uniform temperature above ambient. For 
example, all three curves in Fig. 4(b) nearly overlap and there is barely a noticeable difference (< 0.03°C) in their 
temperature. Also the Ar for the entire SiC remains below the zero difference line between 7:00-15:00 hrs. Even 
though the mirror is warmer than the ambient air by about 2 °C at the beginning, it does not accumulate enough energy 
to raise the temperature of the entire block beyond ambient air until about 15:00 hrs. Under such conditions a mild 
heating of the front surface of SiC mirror during the early morning hours is necessary to bring its temperature close 
to the ambient. 
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2.3. Radial and axial temperature gradients 

Figure (5) shows the day-time and night-time temperature difference of a radial line (r = 0.2 1 m to r = 1 m) along the 
front surface of the Zerodur and SiC mirrors. Temperature evolution is recorded for every two hours interval. Unlike 
SiC, the shape of the temperature curves shown in Fig. 5(a) and Fig. 5(b) for Zerodur mirror is highly distinct for the 
heating and cooling cycles. The gradual onset of solar flux in the morning hours and slow and complex heat diffusion 
across the substrate makes the temperature gradients highly nonuniform. The temperature gradients are particularly 
more pronounced closer to mirror edges. Further, as seen from Fig. 5(b), the predominant heat loss from sides during 
the nigh makes the temperature gradients even more steeper and smoother towards to mirror edges. The numerical 
calculations for the maximum temperature difference (ATmax) existing between the mirror surface and ambient air at 
different times are recoded in the Table [3] 
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Figure 5: Radial temperature difference of the front mirror surface during day-time and night-time for (a)-(b) Zerodur mirror and (c)-(d) for SiC 
mirror. The ambient temperature range is - 10°C. 



Mirror cooling by natural convection is extremely inefficient. So regardless of the material, the temperature 
difference between the mirror surface and surroundings tend to be predominantly high at night, making external 
cooling absolutely necessary. 

To investigate the temperature gradients along the mirror thickness, we chose the axial line (z = 0m to z — 0.203m) 
connecting the three reference points shown in Fig. 1 . Again, the presence of strong and nonuniform axial gradients, 
particularly around the middle of the day and towards the midnight, is clearly evident for Zerodur in Fig. 6(a-b). For 
SiC the net temperature difference varies significantly. However, the apparent lack of axial and radial temperature 
gradients can be seen clearly from Fig. 5(c-d) and Fig. 6 (c-d). 

The radial temperature excess ST-„j and axial temperature excess 6T^ V which we define as the difference in 
maximum and minimum temperature occurring along the radial and axial directions at a particular instant of time, 
is computed for the Zerodur and SiC mirror at different ambient temperature range. These results shown in Table [4] 
indicate that the axial gradients are more severe during the day-time while the radial gradient along the front surface 
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Table 3: Excess mirror surface heating for Zerodur(SiC) mirror under different ambient temperature conditions. 



Time 


Ambient Temperature Range 


(hrs.) 


0°-5°C 


0°-10°C 


0°-20°C 


-10°-0°C 




at f°r i \ 


at f°r i \ 


AT fori 


AT f°(~ i \ 


Uo.UU 


1 .39(-0.17) 


1.56(-0.64) 


0.83(-1.81) 


-0.74(-2.42) 


10:00 


2.16(-0.50) 


0.88(-2.06) 


-1.67(-5.58) 


-0.31(-3.46) 


12:00 


3.03(0.04) 


1.26(-2.11) 


-1.70(-6.61) 


0.54(-3.04) 


14:00 


3.65(0.89) 


1.83(-1.01) 


-0.59(-5.00) 


1.43(-1.68) 


16:00 


3.91(1.82) 


2.72(0.99) 


1.15(-1.26) 


2.17(0.33) 


18:00 


3.42(2.63) 


3.48(3.14) 


3.71(3.60) 


2.81(2.62) 


20:00 


1.73(1.22) 


4.89(4.75) 


7.06(7.62) 


4.29(4.39) 


22:00 


2.55(1.89) 


5.35(4.80) 


8.38(8.28) 


4.84(4.54) 


00:00 


2.75(1.95) 


5.17(4.23) 


8.44(7.51) 


4.79(4.04) 


02:00 


2.61(1.72) 


4.64(3.51) 


7.78(6.33) 


4.36(3.38) 


04:00 


2.31(1.41) 


4.00(2.76) 


6.81(5.04) 


3.80(2.67) 


06:00 


1.98(1.11) 


3.35(2.10) 


5.76(3.86) 


3.21(2.03) 



tend to dominate at night. A general conclusion which follows from a relative comparison of data in Tableland Table 
[4]is that day-time heating leads to 'worse seeing' effects at night. Therefore, a solar telescope mirror cannot be readily 
adapted to carry out the intended astronomical observations at night unless the temperature of the mirror is actively 
controlled. 



Table 4: Numerically computed radial (<5 r rac j) temperature excess and axial temperature excess (ST ax ^) for the Zerodur(SiC) mirror at different 
ambient temperature range. 



Time 


Ambient Temperature Range 


(hrs) 


0" 


-5"C 


0" 


-lO^C 


0"- 


20 U C 


-10" 


-0"C 












*Wc) 






5T^(°C) 


08:00 


0.44(0.01) 


1.64(0.04) 


0.86(0.01) 


1.74(0.04) 


1.26(0.03) 


1.84(0.02) 


0.62(0.04) 


1.72(0.04) 


10:00 


0.42(0.02) 


3.34(0.07) 


0.31(0.04) 


3.42(0.07) 


1.27(0.11) 


3.71(0.06) 


1.23(0.06) 


3.60(0.07) 


12:00 


0.53(0.02) 


3.58(0.08) 


0.48(0.04) 


4.30(0.08) 


2.37(0.14) 


4.74(0.08) 


1.51(0.06) 


4.58(0.08) 


14:00 


0.71(0.02) 


4.15(0.07) 


0.41(0.02) 


4.05(0.07) 


2.61(0.10) 


4.44(0.06) 


1.29(0.03) 


4.34(0.07) 


16:00 


0.94(0.04) 


3.19(0.05) 


0.42(0.02) 


2.92(0.05) 


1.72(0.03) 


2.86(0.04) 


0.76(0.01) 


2.96(0.04) 


18:00 


1.23(0.05) 


1.57(0.01) 


1.05(0.06) 


1.48(0.02) 


0.78(0.07) 


1.22(0.02) 


0.66(0.05) 


1.25(0.01) 


20:00 


1.54(0.06) 


1.15(0.02) 


1.81(0.09) 


1.47(0.02) 


2.11(0.15) 


1.97(0.04) 


1.41(0.08) 


1.22(0.02) 


22:00 


1.67(0.06) 


1.10(0.01) 


2.23(0.10) 


1.62(0.02) 


3.16(0.17) 


2.48(0.04) 


1.91(0.08) 


1.41(0.02) 


00:00 


1.65(0.05) 


1.03(0.01) 


2.44(0.08) 


1.59(0.02) 


3.65(0.15) 


2.54(0.04) 


2.14(0.07) 


1.42(0.02) 


02:00 


1.51(0.04) 


0.90(0.01) 


2.33(0.07) 


1.43(0.02) 


3.73(0.13) 


2.50(0.03) 


2.10(0.06) 


1.29(0.02) 


04:00 


1.33(0.03) 


0.13(0.01) 


2.10(0.05) 


1.23(0.01) 


3.48(0.10) 


2.11(0.03) 


1.93(0.05) 


1.31(0.01) 


06:00 


1.14(0.02) 


0.63(0.01) 


1.83(0.04) 


1.04(0.01) 


3.08(0.08) 


1.78(0.02) 


1.70(0.04) 


0.96(0.01) 



2.4. Thermal deformation 

The conductive flow of heat energy inside a solid material establishes a temperature field which induces thermal 
stresses producing a volumetric changes in the material body. For an isotropic case, the thermal strain s in an elastic 
solid depends on the temperature T and the stress-free reference temperature T ie f which is given by: 

e = a ■ (T - T ref ) (7) 

The const itutive stra in-displacement relationship for an uncoupled thermoelastic problem can be expressed for normal 
strain as (Sadd, 2004J): 

du dv dw 

£ x = —', By = — ; s, = —— (8) 

dx y dy z dz 



and for shear strain as: 




0.05 0.10 0.15 

Mirror thickness (m) 



0.05 0.10 0.15 

Mirror thickness (m) 




(d) 



Mirorr thickness (m) 



0.05 0.10 0.15 
Mirro thickness (m) 



-20:00 

-22:00 

-0:00 

-2:00 

-4:00 

-6:00 



Figure 6: Axial temperature difference along the thickness of the mirror blank during day and night-time for (a)-(b) Zerodur material, and (c)-(d) 
for SiC. The ambient temperature range is - 10°C. 




Figure 7: A snapshot of the surface distortion of (a) Zerodur mirror (b) SiC mirror for (0 - 10°C) ambient temperature change. The direction of the 
displacement vectors at the mirror surface is denoted by arrows shown in green. 



where u, v and w are the displacement components in x, y and z directions, respectively. The FEM numerical imple- 
mentation of Eq.(7)-Eq.(l 1) in the COMSOL Multiphysics software is based on the principle of virtual work done. 
The departure of reflecting surface of the telescope mirror from its nominal shape is computed by calculating the 
displacement p = Vm 2 + v 2 + w 2 of the surface node positions from the undeformed geometry. The discrete sag data 
ds calculated using d s = w — z' Vm 2 + v 2 is particularly useful for representing the surface irregularities in terms of 
Zernike polynomials jPovle et al. , 2002b . For now, we only present a measure of displacement errors across the mirror 
surface computed in terms of the p nns - J pdA/A and w rms = J wdA/A, where A is the area of the reflecting surface 
of the mirror. The optical aberrations are proportional to the rms amplitude, p rms and w rms . Though it does not tell us 
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directly the the shape of the irregularities. 

Fig. 7 shows a graphical snapshot of the front surface distortions of Zerodur and SiC mirror obtained by solving 
structural FEM model. The 3D thermal profiles were supplied as input to the FEM model when the temperature of 
the mirrors had reached maximum value during the day-time heating. In this case, the Zerodur and SiC mirror has 
reached the maximum surface temperature about 12.2°C and 10.3°C, at about 14:40 hrs and 16:30 firs, respectively. 
The large thermal deformations along the mirror edges are greatly due to the curved shape of the mirror surface. 




Radial distance (m) Radial distance (m) 

Figure 8: Calculated displacement along the top surface of (a) Zerodur (b) SiC mirror. 

As discussed in previous section, there are strong thermal gradients within the Zerodur mirror and the resulting 
surface distortions can also be highly irregular. However, due to its extremely low CTE, the magnitude of the thermal 
deformations remains incomparably low. Also, the value of CTE is negative over the operating temperature range. 
The displacement vectors as seen in the Fig. 7(a), are all pointing outward. So the z-displacement and sag, however 
little they may be, are all positive. Though, the near-zero CTE of zerodur is useful to keep the thermal deformation 
at negligible level, the wavefront errors caused by highly nonuniform surface deformations, cannot be corrected by 
simple focus shift. 

For SiC, the displacement vectors are pointing into the mirror. The sag and the z-displacement is, therefore, clearly 
negative. More importantly, these surface perturbations are not uniform over the mirror surface. A typical line profile 
for the total displacement and z-displacement occurring at the peak mirror temperature during the day for Zerodur 
and SiC are shown in Fig. 8. Even though the temperature distribution inside the SiC mirror is highly uniform, the 
total displacement is increasing radially outward from ~9yU to ~ 15yU in a quadric manner. This is not unexpected since 
the original mirror shape itself is parabolic. For a mirror b lank, the lack of temperature gradients in SiC is one of 
the most distinguishing and a highly desirable feature. An uniform temperature rise in SiC mirror does not produce 
significant higher order aberrations. It would only change its global radius of curvature, causing the defocus error in 
the final optical image. The defocus error can be easily corrected by moving the secondary mirror. A smooth quadric 
variations in the surface displacement shown in Fig. 8(b) would support this conclusion. 

The detailed calculations of p rms and w rms at different temperature range for two mirror materials are summarized 
in Table[5] In the FEA structural model, the reference temperature T m f for mirror blanks is assumed to be 20°C. The 
magnitude of thermal distortions would depends on how much the mirror temperature deviates from the reference 
temperature. The calculations shown in Table [5] are therefore, consistent with expected response of the materials 
under specified temperature range. It is also important to note that thermal distortions would vary with the varying 
temperature conditions which requires active or adaptive optics support for corrections. 

3. Conclusion 

The primary mirror heating caused by the light absorption is a complex and challenging problem. The thermal 
performance of a solar telescope mirror needs to be accurately predicted under real observing conditions. This is 
essential to design an efficient and durable temperature control system devised to mitigate the detrimental effects 
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Table 5: The numerical values of the thermal deformation calculated for the front surface of the Zerodur and SiC mirror at different ambient 
temperature range. 



Temp, range 


p rms (jum) 


w rms (jaa) 




Zerodur 


SiC 


Zerodur 


SiC 


-10-0°C 


0.625 


331.28 


0.506 


-261.76 


0-5°C 


0.312 


15.33 


0.260 


-12.13 


0-10°C 


0.216 


10.95 


0.179 


-8.66 


0-15°C 


0.135 


6.95 


0.112 


-5.50 


0-20°C 


0.072 


0.282 


0.058 


-1.81 



of excessive heating during the day. We have outlined an approach to study the thermal and structural response 
of a primary mirror under varying observing conditions. In the FEM model for a 2m class primary mirror, the 
location dependent solar flux and a simple physics based heating and cooling model of the ambient air temperature 
was incorporated. The spatial and temporal evolution of temperature field inside two well known materials for optical 
telescope, Silicon Carbide (SiC) and Zerodur, were examined using 3D numerical simulations. 

The low thermal conductivity of Zerodur mirror gives rise to strong radial and axial temperature gradients that 
are quite distinct for the day-time heating and night-time cooling. Heat loss by free convection is very slow so 
the mirror retain significant heat during the night. The thermal response of the SiC mirror is significantly different 
from the Zerodur. The temperature within the SiC mirror substrate equilibrates rather quickly due to high thermal 
conductivity. The absence of thermal gradients and the advantage of high thermal conductivity of SiC cannot be 
favorably leveraged without some temperature regulation by external means. Thermal distortions of the mirror were 
analyzed using the structural FEM model. High surface distortion seems to result if the operating temperature of the 
mirror deviates significantly from the nominal temperature of the material. In extremely low CTE materials, the mirror 
seeing ultimately limits the telescope performance. It is not the high temperature alone, but the relative incremental 
change in ambient air temperature and the mirror which contribute most to the 'seeing effect'. Large temperature 
changes in ambient air and slow thermal response of the glass materials invariably results in bad seeing. In order to 
augment the scientific productivity of a solar telescope, there have been some suggestions from the astronomers about 
the possibility of utilizing the same telescope for limited, but useful night-time observations. A site with minimum 
diurnal changes in day-night temperature (e.g. 0-5° C) would ideally suit for such dual purpose observations. We 
have only considered the thermal gradients and surface distortions in solid mirror which would be quite different for 
the lightweighted mirror made of same materials. The most conspicuous fallout of the different cell geometries and 
side walls structure is the appearance of thermal footprint in the surface distortions. This would be the subject matter 
of future work. 
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